Flat rotation curves using scalar-tensor theories 
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We computed flat rotation curves from scalar-tensor theories in their weak field limit. Our model, 
by construction, flts a flat rotation profile for velocities of stars. As a result, the form of the 
scalar field potential and DM distribution in a galaxy are determined. By taking into account the 
constraints for the fundamental parameters of the theory (A, a), it is possible to obtain analytical 
results for the density profiles. For positive and negative values of a, the DM matter profile is as 
cuspy as NFW's. 
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PACS numbers: 04.50+h, 04.25.Nx, 98.10.+Z, 98.62. Gq 



I. INTRODUCTION 

Recently, we have worked on some of the effects that 
general scalar-tensor theories (STT) of gravity yield on 
astrophysical scales [H, H, By takingthe weak field 
limit of STT, three parameters appear [J|: the Newto- 
nian constant at infinity Goo, a Compton length-scale 
(A = h/mc) coming from the effective mass of the La- 
grangian potential, and the strength of the new scalar 
gravitational interaction (a). For point-like masses, the 
new Newtonian potential is of a Yukawa type [1, Q and, 
in general, analytical expressions can be found for spher- 
ical [2,] and axisymmetric systems Q. Using these solu- 
tions we build a galactic model that is consistent with a 
flat rotation velocity profile which resembles the obser- 
vations [6]. In the present short contribution, we outline 
a way to construct such a galactic model using STT. 



II. STT AND CONSTRAINTS 

A typical spiral galaxy consists of a disk, a bulge, and a 
dark matter (DM) halo. A real halo is roughly spherical 
in shape and contains most of the matter, up to 90%, of 
the system. Thus, being the halo the main component, 
our galactic model will then consist of a spherical system 
of unknown DM. Therefore, in order compare with ob- 
servations we assume that test particles -stars and dust- 
follow DM particles. To construct our halo model, we 
take the weak field limit of general STT and consider 
a spherically symmetric fluctuation of the scalar fleld 
around some flxed background field. The STT is then 
characterized by a background value of the scalar field, 
(0), a Compton lengthscale. A, and a strength a For 
point-like sources, the new Newtonian potential is well 
known to be of a Yukawa type @ : 
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If one fixes the background field to be the inverse of the 
Newtonian constant, ((/)) — G]^ , then for r ^ A the new 
Newtonian potential coincides with the standard Newto- 
nian one. Given this, for r ^ A one finds deviations of the 



order of (l + a) to the Newtonian dynamics. This setting 
is, however, very constrained by local (solar system) de- 
viations of the Newtonian force, in order for a to be less 
than 10~^° [5]. Alternatively, one can choose the setting 
{(j)) — G^^(l -f a) and the new potential coincides with 
the Newtonian for r ^ A and deviates by 1/(1 -I- a) for 
r ^ A. If one thinks of a galactic system, physical scales 
are around the tens of kiloparsecs -and so is our typi- 
cal A-, then for distances bigger than this, one expects 
constrictions, e.g. from cluster dynamics or cosmology. 
Several authors [3] have considered these deviations, giv- 
ing a rough estimate within the range of —1 < a < 5. 
For example, the value of a = —0.5 yields an asymptotic 
growing factor of 2 in Gn, whereas the value of a = 3.0 
reduces Gn asymptotically by one fourth. 



III. THE GALACTIC MODEL 

For a general density distribution the equations gov- 
erning the weak energy (Newtonian limit) of STT are 
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where (j) is the scalar field fiuctuation and p is the den- 
sity distribution which contains baryons, DM particles, 
or other types of matter. For definiteness we can think 
of DM, but its profile is unknown yet. This will be deter- 
mined by setting the dynamics of test particles, i.e., by 
demanding the new Newtonian potential to be the one 
that solves the rotation curves of test particles -stars and 
dust. That is, we require 



r — - — = const. , 
dr 



(4) 



where this constant is chosen to fit rotation velocities in 
spirals [f)|; this is the simplest model. We then proceed 
to solve for and its solution is: 
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Substituting this result into the original system, ([2l[3]), 
gives 



P = RDM = 



47rGArr2 87r(l + a) 



(6) 



For our convenience, we define the following densities: 
P* = 47rG«r^ ^nd = g^^^+a) ^, wherc p* is the density 
that the system would have to achieve flat rotation curves 
if system would be treated solely with Newtonian physics, 
and p0 is the contribution of the scalar field, which should 
be obtained by integrating the equation: 
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By comparing ^ with ([7]), it seems natural to identify 
m*\/l + a to convert ([7|) to a type ([3]) equation. 



now for m* and p* . The solution is therefore given by 



a=1.0; Ai=1.0 
Pdm/ P / P(t> 
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FIG. 1: All density profiles as a function of the radius are 
shown in a log- log scale. The plots were made with A = 1.0 
and a = 1.0. pdm (solid line), p* (short-dashed line) and 
(long-dashed line). 
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where A* = Ayl + a and a is a length scale of the spher- 
ical system, e.g. related to the distance at which stars 
possess flat rotation curves; the size of the halo is denoted 
by R. Once ^ is solved, its solution is substituted into 
([6]) to find the DM distribution, pdm- The procedure 
outlined here is straightforward. Thus, using our given 
p* the above expression can analytically be integrated to 
have 
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Thus, the DM profile becomes 
vl f 1 

Pdm 
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Following, we plot in figure [T] the density profiles pdm, 
p* and p^ for A = 1, a = 1. The main contribution to 
Pdm comes from p* , which is an inverse squared function 
of the radius. The DM profile is therefore cuspy near 
the galactic centre, similar to the NFW's obtained from 
simulations @. We have computed the best fits to these 
curves for r ^ a, given by pdm ~ ^rm and p 
On the other hand, the best fits for r 3> a are pdm 
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Ln (Pdm) a=l . 0; Ai=0.1; A2 = 0.5, ^3 = 1.0; 1^=5.0 
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FIG. 2: DM density profiles as a function of the radius are 
shown in log- log scale. The plots were made with a? = 1.0 and 
Ai — 0.1 (sofid line), A2 — 0.5 (short-dashed line), A3 = 1.0 
(long-dashed line), and A4 = 5.0 (dashed-dotted Une). 
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and p^ ~ . In figure[2]we plot the DM profile for 
various A values, resulting in small changes in the slope. 
For negative a values, ([9]) becomes negative, but the DM 
profile does not become shallower. We have computed 
the case A — 1.0, a = —0.5 for r <C a and found pdm 
while for r 3> a we obtained pdm 
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IV. CONCLUSIONS 

We have constructed a spherically galactic halo model 
in which we fit the new gravitational potential to match 
a fiat rotation profile for test particles -stars and dust. 
Once this is done, one computes the scalar field fluc- 
tuation and DM profiles resulting from this potential. 
We have found analytical results for our model, Eqs.Q 
and (jlOp . These equations contain the information of the 
scalar field (A, a), which determine the DM profile. We 
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have plotted our results for a = 1.0 and various As. These 
profiles have an inner region as cuspy as the NFW's. Neg- 
ative values of a do not yield a shallower DM profile nei- 
ther. This result may however not surprise since we are 
not resolving the inner structure of the rotation velocity 
curves. We only considered its flat asymptotic behaviour; 



in a forthcoming paper we will study it. 
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